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We study the interaction between several fields initially in coherent states. The solution allows us 
to explain why coherent states remain coherent states when subject to non-Markovian dissipation. 
We first study the interaction between two fields and show that this is the building block of the 
total interaction. We give a completely algebraic solution of this system. 

I. INTRODUCTION 

66 '. 

It is well known that a coherent state subject to dissipation keeps its form during the dynamics. This is, given the 
' master equation for a field in a lossy cavity at zero temperature 

<n : 

Ch -p = ^(apa) — a)ap — pa) a) (I) 



with a {a)) the annihilation (creation) operator for the cavity mode, 7 the decay constant and p the density matrix, 

' if the initial state of the cavity field is a coherent state |cv), then the dynamics shows that it will decay in time as 

\ae rt ) (see for instance [![. One possible answer about why the coherent states preserve its form during decay is 

the fact that coherent states are eigenstates of the annihilation operator, however this argument does not hold for a 

'"5 ' dissipative two-photon process @ 
ft 
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t = 7(aV 2 -a t2 aV-^ t2 « 2 ) (2) 

a- dt ' 

even though coherent states are also eigenstates of the annihilation operator squared (so do even and odd coherent 
states 0). 

Both equations above are obtained using Born-Markov approximations 0, E[ • In the case in which such approxi- 
mations are not used, i.e. when the interaction between a harmonic oscillator and a set of harmonic oscillators (the 
environment) is considered, is not clear how a coherent state decays. Here we will try to answer this question. First 
we will consider the interaction between two fields, to later generalize the result to a field interacting with many. 



II. TWO FIELDS INTERACTING 



o . 

OO , Consider the Hamiltonian of two interacting fields (we set h = 1) 

o ■ 

H = u a a f a + u b b ] b + A(a f 6 + 6 f a), (3) 

a system like this may be produced by the interaction of two quantized fields with a two-level atom Q by transforming 
?H ' to the interaction picture, i.e. getting rid off the free Hamiltonians, we obtain 

a ; 

= Aa t fl + A(a t 6 + 6 t a), (4) 
with A = u a — uji,, the detuning. It is useful to define normal- mode operators by Q 

Ai = da + 76, A 2 = 7a - Sb, (5) 

with 



A 2A n-A 

S= ^m-Ay 7 = V^r (6) 

with O = V A 2 + 4A 2 the Rabi frequency. A\ and A2 are annihilation operators just like a and b and obey the 
commutation relations 

[a 1j V] = [a 2 ,4] = i, (7) 
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moreover, the normal-mode operators commute with each other 

[A U A 2 ] = [Ai,4] = 0. 
In terms of these operator the Hamiltonian Q becomes 

Hi = fiiA{Ai + ii 2 A\A 2 , 



(8) 



(9) 



with ni : 2 — (A ± f2)/2. Up to here, we have translated the problem of solving Hamiltonian (1) into the problem 
of obtaining the initial states, for the "bare" modes in the initial states for the normal modes. In order to have a 
way of transforming states from one basis to the other, we note that the vacuum states in both systems |0) a |0)(, and 
|0)ai|0)a 2 differ only for a phase [6J. First note that 



Ai|0) a |0) b = 0, 



(10) 



and in a similar way it may be seen the other normal-mode annihilation operator, A 2 , has the same effect. Choosing 
the phase so that 



|0) a |0) 6 = |0) Al |0U 2 . 

If we consider coherent states as initial states for the interaction, we obtain the evolved wavefunction 

|^(t)) = e- u ^ A ^ A ^D a (a)D b (^) a \0) bl 
= e- lt ^ A ^ A ^D a {a)D b {P)\Q) Al |0)a 2 



(11) 



(12) 



where the D c (e) = exp(ec^ — e*c) is the Glauber displacement operators |7]. From (JSJ) we can write the operators o 
and b in terms of the operator A\ and A 2 (fT2"|) as 



W)) = e- lt ^ A ^+^ A ^D Al (aS + P^DA^aj - P5)\0) Ai \0)a 2 - (13) 
Passing the exponential in the above equation to the right and applying it to the vacuum states we obtain 

m)) = D Al ([a5 + H^^DaAWi - P5]e- lll2t mAMA 2 

= {[aS + ^e-^U^-^e-^A,. (14) 

Equation (|14p shows that in the new basis, coherent states remain coherent during evolution. By transforming back 
to the original basis we obtain 



\i/>(t)) = \S[aS + /^e"^ 1 * + 7 [a 7 - (36}e- lt * 2t ) a \~/[a5 + /3 7 ] e - Vlt - S[aj - (3S}e'^ 2t ) b , 



(15) 



i.e. coherent states remain coherent during evolution. This will be used next Section as the building block for the 
interaction of many modes. In obtaining (fT4")) and (fT4"|) , we have used the following property 



D c ( £l )D c (e 2 ) = £> c (ei + £l ) e 3(^-^= 



(16) 



III. GENERALIZATION TO n MODES 



Consider the Hamiltonian of the interaction of k fields 



From the Hamiltonian above, we can produce the following matrix 



/ wi A 2 i 

A12 UJ 2 

M Al3 A23 



\A lfe A 



2k 



Afei \ 

Afc2 
Afc3 



(17) 



(18) 
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i.e. D is a diagonal matrix whose elements are the eigenvalues of the matrix M, defined from the Hamiltonian. The 
matrix R is therefore M's eigenvectors matrix. The solution to the Schrodinger equation subject to the Hamiltonian 
(fTTf with all the modes initially in coherent states, \ip(0)) = \ax)x\ot2)2---\oi n ) n , is simply the direct product of coherent 
states 

m)) = \n ■ P(t))i\?2 • P(t)) 2 ...\r n ■ 0{t)) n (29) 

with f3(t) — {f\ ■ ae _ ' Mli , r 2 • d?e~ V2 *, r n ■ ae - '* 1 "') and the vector a — {ax, ev 2 , a n ) is composed by the coherent 
amplitudes of the initial wave function. Up to here we have shown that the interaction of several modes initially in 
coherent states does not change the form of those states (remain coherent), but modifies their amplitude. If we choose 
the interaction constants to be Xxj 7^ for 1 ^ j and the rest as zero, we are dealing with the interaction between 
one field and n — 1 fields. If n — > 00 and the amplitudes otj are zero for j > 1, we deal with the interaction of one 
field with n — 1 one of them in a coherent state with amplitude ax and the rest in the vacuum. Therefore, the most 
likely situation we have is the coherent state decaying towards the vacuum while keeping its coherent form. 



IV. CONCLUSIONS 



We have shown that a system of n interacting harmonic oscillators initially in coherent states, remain coherent 
during the interaction. In particular, if one considers one field (harmonic oscillator) interacting with many fields 
(harmonic oscillators), i.e. consider only Xxj 7^ and Xjx 7^ for j > 2, and all the others to be zero, we can model 
non-Markovian system-reservoir interaction. If we consider the system to be in a coherent state and all the others 
fields that form the environment in a vacuum state (this is also in coherent states with zero amplitude), after evolution, 
the amplitude of the coherent state will diminish, as one photon will go to another mode, keeping its coherent nature. 
If the number of modes that form the environment is very large, an event of the photon going back to the system is 
quite unlikely. Therefore the next probable event is precisely the loss of another photon by the system, etc. until it 
arrives to a state close to the vacuum. In case the number of modes interacting with the system is infinite, then the 
vacuum would be the final state of the system. In other words, the total system perform the following transition 

|a)i|0) 2 ...|0) n ^|<5 1 ) 1 |5 2 ) 2 ...|5 2 ) n , (30) 

where the coherent amplitudes,^ — > 0, as n — > 00. 

In conclusion we have given a complete algebraic solution to the problem of n interacting harmonic oscillators, 
without Born-Markov approximations. 
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